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ABSTRACT 



We estimate the mass density fluctuations power spectrum (PS) on large scales 
by applying a maximum likelihood technique to the peculiar velocity data of the 
recently completed redshift-distance survey of early- type galaxies (hereafter ENEAR) . 
■ Parametric CDM-like models for the PS are assumed, and the best fit parameters are 

determined by maximizing the probability of the model given the measured peculiar 
• velocities of the galaxies, their distances and estimated errors. The method has been 

applied to CDM models with and without COBE normalization. The general results 
are in agreement with the high amplitude power spectra found from similar analysis 
of other independent all-sky catalogs of peculiar velocity data such as MARK III 
(Willick et al. 1997)and SFI (Giovanclli et al. 1998; da Costa et al. 1996), in spite of 
the differences in the way these samples were selected, the fact that they probe different 
regions of space and galaxy distances are computed using different distance relations. 
For example, at k = 0.1 /iMpc -1 the power spectrum value is P(k)fl 1 ' 2 = (6.5 ± 3) x 
10 3 ( /i -1 Mpc) 3 and rjs = er^ ' 6 = 1-1^0 35' the quoted uncertainties refer to 3cr error 
level. We also find that, for ACDM and OCDM COBE-normalized models, the best-fit 
parameters are confined by a contour approximately defined by Qh 13 = 0.377 ± 0.08 
' and fi/i 88 = 0.517 ±0.083 respectively. T-shape models, free of COBE normalization, 

results in the weak constraint of T > 0.17 and in the rather stringent constraint of 
778 = 1.0 ± 0.25. All quoted uncertainties refer to 3cr confidence-level (c.L). 

The calculated PS has been used as a prior for Wiener reconstruction of the density 
field at different resolutions and the three-dimensional velocity field within a volume of 
radius « 80/i _1 Mpc. All major structures in the nearby universe are recovered and are 
well matched to those predicted from all-sky redshift surveys. The robustness of these 
features has been tested with Constrained Realizations (CR). Analysis of reconstructed 
three-dimensional velocity field yields a small bulk flow amplitude (~ 160 ±60 kms -1 
at 60 ft^Mpc) and a very small rms value of the tidal field (~ 60 kins -1 ). The results 
give further support to the picture that most of the motion of the Local Group arises 
from mass fluctuations within the volume considered. 

Key words: cosmology: observations - cosmology: theory - dark matter - galaxies: 
distances and redshifts - large-scale structure of universe - methods: statistical 



1 INTRODUCTION 

The canonical model of cosmology assumes that large-scale 
structure has grown out of small density perturbations via 
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the process of gravitational instability. These initial fluc- 
tuations are usually assumed to satisfy the statistics of a 
Gaussian random field, solely characterized by its power 
spectrum. In the linear regime, the fluctuations grow self- 
similarly and retain their initial distribution and power spec- 
trum shape. Therefore, mapping the underlying cosmolog- 
ical velocity field and its power spectrum on large scales, 
provides a direct probe to the origin of structure in the uni- 
verse. 

The PS, the three-dimensional distribution of luminous 
matter and the predicted peculiar velocity field have been 
derived from a variety of data sets, especially from all-sky 
redshift surveys (e.g. Strauss & Willick 1995 for a review of 
earlier work; Sutherland et al. 1999; Branchini et al. 1999). 
Unfortunately, however, the distribution of galaxies in these 
catalogs is not necessarily an unbiased tracer of the underly- 
ing mass distribution, and suffer from the infamous "galaxy 
biasing" problem. Furthermore, in estimates from redshift 
surveys, uncertainties arise from the complicated relation 
between the real space and the redshift space distributions, 
known as redshift distortions (e.g., Kaiser 1987, Zaroubi and 
Hoffman 1996). In order to avoid these problems altogether 
it is advantageous to appeal to dynamical data, in particular 
catalogs of galaxy peculiar velocities on large scales. 

Peculiar velocities enable a direct and reliable deter- 
mination of the mass PS and distribution, under the nat- 
ural assumption that the galaxies are unbiased tracers of 
the large-scale, gravitationally-induced, velocity field. Fur- 
thermore, since peculiar velocities are non-local and have 
contributions from different scales, analysis of the peculiar 
velocity field provides information on scales somewhat larger 
than the sampled region (e.g. Hoffman et al. 2000). For the 
same reason peculiar velocities are adequately described by 
linear theory even when densities become quasi-linear (e.g., 
Freudling et al. 1999). Consequently, the dynamics and the 
distribution of peculiar velocities are well described by the 
linear regime of gravitational instability and by a Gaussian 
probability distribution function (PDF), respectively. 

Assuming that both the underlying velocity field and 
the errors are drawn from independent random Gaussian 
fields, the observed peculiar velocities constitute a multi- 
variant Gaussian data set, albeit the sparse and inhomo- 
geneous sampling. The corresponding posterior PDF is a 
multivariate Gaussian that is completely determined by the 
assumed PS and covariance matrix of errors. Under these 
conditions one can write the joint PDF of the model PS and 
the underlying velocity or density field. 

The purpose of the present study is to calculate, from 
the joint PDF, the PS and 3D mass distribution, as well as 
the 3D peculiar velocity field, as derived from the newly com- 
pleted ENEAR galaxy peculiar velocity catalog (da Costa 
et al. 2000a, Paper I). First, the PS model parameters are 
estimated by maximizing the likelihood function given the 
model (Zaroubi et al. 1997). An identical likelihood estima- 
tion of the power spectrum has been previously applied to 
the Mark III (Zaroubi et al. 1997) and the SFI (Freudling 
et al. 1999) data sets. In both cases the analysis yielded a 
high amplitude power spectrum. Although the results from 
those two catalogs are consistent with each other, they are 
marginally inconsistent with the power spectra measured 
from redshift catalogs (e.g., da Costa et al. 1996; Suther- 
land et al. 1999), inferred from the analysis of the velocity 



correlation function (e.g. Borgani et al. 2000a, 2000b), and 
from velocity-velocity comparisons (e.g., Davis et al. 1996, 
da Costa et al. 1998). One of our goals is to use the same 
methodology employed before for the Mark III and SFI to 
the new ENEAR catalog to directly test the reproducibility 
of the results with an independent sample based on a differ- 
ent distance indicator but probing a comparable volume. 

Second, the Wiener filter (WF) solution of the field is 
recovered by finding the most probable field given the PS 
and the data (Zaroubi et al. 1995, 1999). Constrained real- 
izations (CR) are then used to sample the statistical scatter 
around the WF field (Hoffman and Ribak 1991). The mass 
density PS is used to calculate the smoothed Wiener filtered 
density and 3D velocity fields given the measured radial ve- 
locities (Zaroubi et al. 1995, 1999). The WF provides an 
optimal estimator of the underlying field in the sense of a 
minimum-variance solution given the data and an assumed 
prior model (Wiener 1949; Press et al. 1992). The prior 
defines the data auto-correlation and the data-field cross- 
correlation matrices. In the case where the data is drawn 
from a random Gaussian field, the WF estimator coincides 
with the conditional mean field and with the most proba- 
ble configuration given the data (see Zaroubi et al. 1995). 
It should be noted that Kaiser & Stebbins (1991) were the 
first to propose a Bayesian solution to the problem of recon- 
struction from peculiar velocity data sets. 

Finally, the recovered three-dimensional velocity field is 
used to compute the amplitude of the bulk flow and to de- 
compose the velocity field in terms of a divergent and tidal 
components which enables one to separate the contribution 
to the measured peculiar velocity field from mass fluctu- 
ations within and outside the volume probed by the data 
(Hoffman et al. 2000). 

The methods adopted in this study do not involve any 
explicit window function, weighting or smoothing the data. 
In addition, they automatically underweight noisy, unreli- 
able data. However, a few simplifying assumptions are re- 
quired: 1) peculiar velocities are drawn from a Gaussian ran- 
dom field; 2) peculiar velocities are related to the densities 
through linear theory; 3) errors in the D n — a inferred dis- 
tances constitute a Gaussian random field with two compo- 
nents, the first scales linearly with distance while the second 
models the nonlinear evolution of the velocities as a constant 
scatter. The need to assume a parametric functional form for 
the PS is also a limitation. 

The outline of this paper is as follows. In § ^ we briefly 
describe the peculiar velocity data used in the present anal- 
ysis. The PS analysis is carried out in § a. The Wiener filter- 
ing is applied to the ENEAR data in § H where maps of the 
density field are presented and compared to those predicted 
from redshift surveys. Also shown in this section are the re- 
covered three-dimensional velocity field and the results of its 
analysis. Our results are summarized and discussed in § El 



2 THE DATA 

In the present analysis, we use the ENEAR redshift-distance 
survey described in greater detail in Paper I of this se- 
ries. Briefly, the ENEAR sample consists of roughly 1600 
early- type galaxies brighter than tub = 14.5 and with 
cz < 7000 kms -1 , for which D n — odistances are avail- 
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able for 1359 galaxies. Of these 1145 were deemed suitable 
for peculiar velocity analysis according to well-defined cri- 
teria (Paper I; M. Alonso et al. , in preparation). To the 
magnitude- limited sample we added 285 fainter and/or with 
redshifts > 7000 km s -1 , 129 within the same volume as the 
magnitude-limited sample. In total, the cluster sample con- 
sists of 569 galaxies in 28 clusters, which are used to derive 
the distance relation (M. Bernardi et al. , in preparation). 
Over 80% of the galaxies in the magnitude-limited sample 
and roughly 60% of the cluster galaxies have new spectro- 
scopic and R-band photometric data obtained as part of 
this program. Furthermore, repeated observations of several 
galaxies in the sample (M. Alonso et al. , in preparation; 
G. Wegner et al. 2000) provide overlaps between observa- 
tions conducted with different telescope/instrument config- 
urations and with data available from other authors. These 
overlaps are used to tie all measurements into a common sys- 
tem, thereby ensuring the homogeneity of the entire dataset. 
In contrast to other samples new observations conducted by 
the same group are available over the entire sky. The com- 
parison between the sample of galaxies with distances and 
the parent catalog also shows that the sampling across the 
sky is uniform. 

Individual galaxy distances were estimated from a di- 
rect D n — a template relation derived by combining all the 
available cluster data (M. Bernardi et al. , in preparation), 
corrected for incompleteness and associated diameter-bias 
(Lynden-Bell et al. 1988). From the observed scatter of the 
template relation the estimated fractional error in the in- 
ferred distance of a galaxy is A ~ 0.19, nearly independent 
of the velocity dispersion. 

Since early-type galaxies are found preferentially 
in high-density regions, galaxies have been assigned to 
groups/clusters using well-defined criteria imposed on their 
projected separation and velocity difference relative to the 
center of groups and clusters. These systems were iden- 
tified using objective algorithms applied to the available 
magnitude-limited samples, comprising all morphological 
types, with complete redshift information probing the same 
volume. For membership assignment we used group cata- 
logs published by Geller & Huchra (1983), Maia, da Costa 
& Latham (1988), Ramella, Pisani & Geller (1997) as well as 
unpublished results (M. Ramella et al. , in preparation) cov- 
ering other regions of the sky. The characteristic size and ve- 
locity dispersion of these groups/clusters were used to estab- 
lish the membership of the ENEAR early-types, as described 
in Paper I. We find isolated galaxies, groups with only one 
early-type, and groups with two or more early-types. Early- 
type galaxies in a group/cluster are replaced by a single 
object having: (1) the redshift given by the group's mean 
redshift, which is determined considering all morphologies; 
(2) the distance given by the error-weighted mean of the in- 
ferred distances, for groups with two or more early-types; 
and (3) the fractional distance error given by A/y^jV), 
where N is the number of early-types in the group. In some 
cases groups were identified with Abell/ACO clusters within 
the same volume as the ENEAR sample and fainter cluster 
galaxies were added, as described in Paper I. In the analysis 
below we compute the dipole component of the velocity field 
out to 6000 kms~ as probed by all objects, and by splitting 
the sample into two independent sub-samples consisting of 
field galaxies and groups/clusters. The latter is done to eval- 



uate directly from data the amplitude of possible sampling 
errors. 

The inferred distances are corrected for the homoge- 
neous and inhomogeneous Malmquist bias (1MB). The lat- 
ter was estimated using the PCSz density field (Branchini et 
al. 1999), corrected for the effects of peculiar velocities, in 
the expressions given by Willick et al. (1997). In this cacu- 
lation we also include the correction for the redshift limit of 
the sample. A complete account of the sample used and the 
corrections applied will be presented in a subsequent paper 
of this series (M. Alonso et al. , in preparation). 



3 POWER SPECTRUM 

The calculation of the matter PS from the peculiar velocity 
data by means of likelihood analysis requires a relation be- 
tween the velocity correlation function and the power spec- 
trum. Define the two-point velocity correlation (3 x 3) tensor 
by the average over all pairs of points Vi and Tj that are sep- 
arated by r = Tj — Ti , 



4V(r) = (*v(ri)v„(rj)), 



(1) 



where « M ( r i s the n component of the peculiar velocity 
at Ti. In linear theory, it can be expressed in terms of two 
scalar functions of r = |r| (Gorski 1988), computed from 
the parallel and perpendicular components of the peculiar 
velocity, relative to the separation vector r, 



1V(r) = *x(r)5 M „ + [*||(r) - *x(r)]f M rv. 



(2) 



The spectral representation of these radial correlation func 
tions is 



2tt 2 



Jo 



P{k)K ±;ll {kr) dk, 



(3) 



where K±(x) = ji(x)/x and K\\(x) = jo — 2ji(x)/x, with 
ji{x) the spherical Bessel function of order /. The cosmo- 
logical Q dependence enters as usual in linear theory via 
/(f2) « O ' 6 , and Ho is the Hubble constant. A paramet- 
ric functional form of P(k) thus translates to a parametric 
form of $,,1,. Note that the quantity that can be derived 
from peculiar- velocity data via the linear approximation is 
/ 2 (fi)P(fc), where P is the mass density PS. 

Let m be the vector of model parameters and d the 
vector of A?" data points. Then Bayes' theorem states that 
the posterior probability density of a model given the data 



7>(m|d) 



P(m)P(d|m) 
P(d) ■ 



(4) 



The denominator is merely a normalization constant. 
The probability density of the model parameters, V(m), is 
unknown, and in the absence of any other information we 
assume it is uniform within a certain range. The conditional 
probability of the data given the model, - P(d|m), is the like- 
lihood function, £(d|m). The objective in this approach, 
which is to find the set of parameters that maximizes the 
probability of the model given the data, is thus equivalent 
to maximizing the likelihood of the data given the model 
(cf. Zaroubi et al. 1997; Jaffe & Kaiser 1994). The Bayesian 
analysis measures the relative likelihood of different mod- 
els. An absolute frequentist's measure of goodness of fit is 
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provided by the x 2 P er degree of freedom (hereafter, d.o.f.), 
which we use as a check of the best parameters obtained by 
the likelihood analysis. 

Assuming that the velocities are a Gaussian random 
field, the two-point velocity correlation tensor ^ fully char- 
acterizes the statistics of the velocity field. Define the radial- 
velocity correlation (N x N) matrix Uy by 



Uij = fif^fj = ±(r)sin8isin8j + 'Py (r)cos#iCos#j 



(5) 



where i and j refer to the data points, r = |r| = \rj — r; 
and the angles are defined by 8i = rVr (Gorski 1988; Groth, 
Juszkiewicz & Ostriker 1989). Let the inferred radial pecu- 
liar velocity at ri be u°, with the corresponding error ti also 
assumed to be a Gaussian random variable. The observed 



correlation matrix is then Uij 
hood of the TV data points is 



Ui, 



e t Sij , and the likeli- 



£ = [(2^)^ det(f/y)] 



-1/2 



exp 



(6) 



and positive definite, we can use the Cholesky decomposi- 
tion method (Press et al. 1992) for computing the likeli- 
hood function (Eq. (^i|)). The significant contribution of the 
errors to the diagonal terms makes the matrix especially 
well-conditioned for decomposition. 

The errors are assumed to have two contributions, the 
first is the usual D n — a distance proportional errors (about 
19% per galaxy for ENEAR). The second is a constant er- 
ror that accounts for the non-linear velocities of galaxies in 
the high density environment in which early-type galaxies 
reside. This term represents our poor understanding of the 
complex correlations introduced by non-linear evolution. For 
each power spectrum model, we have performed the likeli- 
hood analysis assuming this constant value to be either null 
or 250 kms -1 but, as shown below, the difference in the 
results are only marginal and do not affect our general con- 
clusions. 



3.1 COBE-Normalized CDM Models 

Wc first restrict our attention to the generalized family of 
CDM cosmological models, allowing variations in the cos- 
mological parameters £1, A and h. Furthermore, four-year 
COBE normalization is imposed as an additional external 
constraint. The general form of the PS for these models is 



P(k) = A CO BE(n, Q, A) T 2 (n, Q B , h; k) k" 



(7) 



where the CDM transfer function proposed by Sugiyama 
(1995) is adopted, 



T(k) 



k 



ln(l + 2.3g)) 
2.34g 

[l + 3.89g + (16.1g) 2 + (5.46g) 3 + (6.71g) 4 
fi/l exp(-0(, - hlJa^lb/VL) (hMpc" 



C 



-1/4 



(9) 



The parameters Q and h are varied such that they span the 
range of currently popular CDM models, including ACDM 
(fi + A = 1, n < 1) and OCDM (A = 0, Q < 1). In all cases, 
the baryonic density is assumed to be Qt = 0.019ft -2 , which 
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Figure 1. Contour map of ln-likclihhod in the h — Q plane for 
ACDM models with 250 kms -1 thermal error component (upper 
panel) and with zero thermal error (lower panel). The contours 
denote the most likely values within 1,2 and 3<r c.l.. 



is the value currently favored by primordial nucleosynthesis 
analysis (e.g., Buries & Tytler 1998). We limit our inquiry to 
models without tilt, namely to models with n = 1. For each 
model, the normalization of the PS is fixed by the COBE 4- 
year data (Bennet et al. 1996); for more details see Zaroubi 
et al. (1997 and references therein). 

Figure la shows the likelihood contour map in the Q — h 
plane, for the ACDM family of models with n = 1 (nor- 
malization by Sugiyama 1995), the error matrix is assumed 
here to have an additional diagonal random contribution of 
250 km s" 1 that accounts for the nonlinear evolution of the 
galaxies. The most probable parameters in this case (in the 
range Q, < 1) are SI = 1 and h = 0.5. The elongated contours 
clearly indicate that neither Q nor h are independently well 
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Figure 2. Same as in Figure la but for OCDM models. 



constrained. It is rather a degenerate combination of the 
two parameters, approximately £lh x with x ~ 1 (i.e. a com- 
bination close to the V parameter) that is being determined 
tightly by the elongated ridge of high likelihood. 

Figure lb shows the likelihood results for the same 
ACDM model shown in Fig. la but with no random con- 
tribution to the error matrix. The contours in Fig. lb show 
very little changes relative to those shown in panel (a), no- 
tably they get tighter and the best values of Q for a given 
Hubble constant are somewhat higher. The addition of a 
reasonable random component to the error matrix does not 
alter the results in any significant way for any of the PS 
models considered in this study. For the rest of the PS mod- 
els we show the calculation with the addition of a constant 
error of 250 km s -1 . 

Figure 2 shows the similar likelihood map for OCDM 
with n = 1. The most probable values here are f2 = 0.53 
and h — 1. The values of Q and h are not independently 
constrained here as well. 

We can thus quote stringent constraints on the con- 
ditional best value of fi given h for the COBE normal- 
ized CDM models shown in Figs, la and 2: Q ~ (0.377 ± 
0.08)/i" 13 for ACDM, and Q « (0.517 ± 0.083)/i" 88 for 
OCDM. 



3.2 The T Model 

To recover the PS from the velocity data independent of 
the COBE normalization, we use as a parametric prior the 
so-called F model (e.g., Efstathiou, Bond and White 1992), 



P(k) = AkT 2 (k), 

T(k) = (i + K/r + (Wc/r) 3/2 -Kcfc/rfrJ (10) 

with a = 6.4/i _1 Mpc, b = 3.0/i _1 Mpc, c = 1.7/i _1 Mpc 
and v = 1.13. The free parameters to be determined by the 
likelihood analysis are the normalization factor r]s = asQ. ' 6 
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Figure 3. Contour map of ln-likelihood for the T model in the 
r — r;8 plane. The contours denote the 1, 2 and 3<r 

c.l. 

and the V parameter. In the context of the CDM cosmo- 
logical model, F has a specific cosmological interpretation, 
T = flh. Here, however, Eq. ( Jio) ) serves as a generic func- 
tion with logarithmic slopes n = 1 and —3 on large and 
small scales respectively, and with a turnover at some inter- 
mediate wavenumber that is determined by the single shape 
parameter V. 

Figure 3 shows the contour map of ln£ in the P-fys 
plane. Although the likelihood analysis poses strong con- 
straint on the allowed values of r]s (= 28 w ith 3<r c.l.), 
it only weakly constrains the value of F (> 0.18 with 3<r 
c.l.), and r = 0.25 is excluded with 2<r c.l.. 



3.3 Results and Comparison between the Various 
Models 

The best fit models for each CDM family have a comparable 
likelihood, with the most likely model is the OCDM model 
with Q — 0.53 and h = 1. All best fit models agree within 
w 20% for k > 0.1/iMpc~\ The amplitude of the PS at 
k = 0.1 ftMpc" 1 for all models lies within P(k)n 12 = (6.5± 
3) x 10 3 ( /i _1 Mpc) 3 and the values of r]8 are within the range 
1 l +0 ' 2 

Figure 4 shows the power spectrum of the most-likely 
COBE-normalized model and the 3cr errors about it. It also 
shows the PS corresponding to most likely models of ACDM 
and r models. Within the errors, the most likely PS for each 
CDM family are very consistent, especially at intermediate 
scales (30 — 50/i~ 1 Mpc), where the data information con- 
tent provide the strongest constraint. Also shown in Fig. 4 
are the best-fit PS, obtained from a similar likelihood anal- 
ysis, for the Mark III and SFI data sets. As can be seen 
the most likely PS for the three catalogs are in good agree- 
ment. This result shows that the high amplitude PS found 
from peculiar velocity data is unlikely to be due to possible 
non-uniformities of these catalogs or to the type of galaxies 
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Figure 4. The PS of the most probable COBE-normalized 
OCDM (solid-bold) and ACDM (dashed-bold) models and of 
the T-model (dotted-dashed-bold). shown also the most proba- 
ble models as estimated from Mark III (dotted-dashed) and SFI 
(triple-dotted-dashed) data sets. The shaded region around the 
PS marks the 3-cr c.l. note that the dynamical range of the data 
is confined to 0.05 k ^5 0.3. 



used. In fact, while Mark III and SFI relied predominantly 
on TF distances to spirals, ENEAR relies on D n — a dis- 
tances to early-type galaxies. On the other hand, the rea- 
son for the discrepancy in the cosmological constraints be- 
tween the maximum likelihood method and other methods 
(da Costa et al. 1998; Strauss & Willick 1998; Borgani et al. 
2000a, 2000b) remains unresolved. The former yields a sys- 
tematically higher amplitude PS, as reflected by the high 
values of rjs, which is also in disagreement with the con- 
straints derived from other analysis of LSS data. Possible 
explanations are given in § [s| 

In all the COBE-normalized PS models considered the 
X 2 /d.o.f. of the best fit models is of the order of 0.93. This 
value deviates by about 2a from the \ 2 /d.o.f desired value 
of unity. This, however, does not pose any serious problem 
since many of the models within the likelihood most likely 
contours have a x 2 /d.o.f . ~ 1. The \ 2 /d.o.f for the r- model 
is 0.99. 



4 WIENER FILTER & CONSTRAINED 
REALIZATIONS 

4.1 The Method 

Having determined the power spectrum, all the ingredient 
needed to Wiener reconstruct the density and velocity fields 
are ready. Details on the general application of the WF/CR 
method to the reconstruction of large-scale structure are de- 
scribed in Zaroubi et al. (1995), where the theoretical foun- 
dation is discussed in relation with other methods of estima- 
tion, such as Maximum Entropy. The specific application of 
the WF/CR method to peculiar velocity data sets has been 



presented in Zaroubi et al. (1999). Here we provide only a 
brief description of the WF/CR method, for more details the 
reader is referred to the original references references given 
above. 

We assume that the peculiar velocity field v(r) and 
the density fluctuation field 5(r) are related via linear 
gravitational-instability theory. Under the assumption of a 
specific theoretical prior for the power spectrum P(k) of the 
underlying density field, one can write the WF minimum- 
variance estimator of the fields as 



WF/ \ 

v (r) 



and 



Ui Uj 



* WF (r) = (J(rK )(«?«; 



(11) 



(12) 



A well known problem of the WF is that it attenuates 
the estimator to zero in regions where the noise dominates. 
The reconstructed mean field is thus statistically inhomo- 
geneous. In order to recover statistical homogeneity we pro- 
duce constrained realizations (CR), in which random realiza- 
tions of the residual from the mean are generated such that 
they are statistically consistent both with the data and the 
prior model (Hoffman and Ribak 1991; see also Bertschinger 
1987). In regions dominated by good quality data, the CRs 
are dominated by the data, while in the limit of no data the 
realizations are practically unconstrained. 

The CR method is based on creating random realiza- 
tions, <5(r) and v(r), of the underlying fields that obey the 
assumed PS and linear theory, and a proper set of random 
errors e, . The velocity random realization is then "observed" 
like the actual data to yield a mock velocity dataset u°. 
Constrained realizations of the dynamical fields are then ob- 
tained by 



v CR (r) = v(r) + ( v(r)w° 



and 



5 CK (r) = S(r) + /5(r)u°\/u°u°\ (u° - u°) 



(13) 



(14) 



The two types of covariance matrices in the above equations 
are computed within the framework of linear theory as fol- 
lows. The covariance matrix of the data < u°u° > is the same 



matrix Uij that appears in eq. |6J. 

The cross-correlation matrix of the data and the under- 
lying field enters the above equations as, e.g., 



5(t)u°) = («(r)vfo)) 



(15) 



The two-point cross-correlation vector between the density 
and velocity fields is related to the PS via 



<5(x)v(x + r] 



Homo) 

2tt 2 



kP(k)j 1 (kr)dk. (16) 



The assumption that linear theory is valid on all scales 
enables us to choose the resolution as well, and in particular 
to use different smoothing radii for the data and for the re- 
covered fields. In our case no smoothing were applied to the 
radial velocity data while we choose to reconstruct the den- 
sity field with a finite Gaussian smoothing of radius R. This 
alters the density-velocity correlation function by inserting 
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SGX (Mpc/h) 

Figure 5. Reconstruction from the ENEAR catalog. Left panel: 
The density field map in the Supergalactic plane, with G12 
smoothing. Density contour spacing is 0.1, positive contours are 
solid, negative contours are dashed and 8 = is denoted by heavy- 
solid line. The shading indicates regions where the error is less 
than 0.3. 

the multiplicative term exp[— k 2 R 2 /2] into the integrand of 
eq. 

A theoretical estimate of the signal-to-noise ratio (S/N) 
at every point in space is given by a simple expression (see 
Zaroubi et al. 1999) but it requires the calculation and inver- 
sion of very large matrices. Therefore, in this study we esti- 
mate the point to point error by conducting a large number 
of CRs. In the case of random Gaussian fields, the ensemble 
of CRs defined in eq. (|l^)and eq. (|l4|) samples the distribu- 
tion of uncertainties in the mean Wiener density and velocity 
fields (Hoffman & Ribak 1991). 

It is worth noting that the WF represents a general 
minimum-variance solution under the sole assumption that 
the field is a random field with a known power spectrum. 
No assumption has to be made here regarding higher order 
correlations (or the full joint probability distribution func- 
tions) of the underlying field. On the other hand, the CRs 
are derived under the explicit assumption of a full Gaussian 
random field. 

4.2 Maps of Density and Velocity Fields 

Figure 5 shows the map of the density field along the Super- 
galactic plane obtained from the ENEAR data using a Gaus- 
sian smoothing radius of 1200 kms -1 (hereafter G12). The 
shaded area corresponds to the region where the error, as es- 
timated from performing 10 CRs, in density is less than 0.3. 
The main features of our local universe are easily identified 
in the WF map, including the Great Attractor (GA) on the 
left and the Perseus-Pisces supercluster (PP) in the lower 
right. There is also a hint of the Coma cluster, which lies 
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Figure 6. Upper panel: T^ samePas in Fig. 5 but with G9 
smoothing. The shaded area indicates regions with error smaller 
than 0.45. Lower panel: The G9 smoothing density reconstruction 
from the PSCz redhsift catalog (Branchini et al. 1999) 



just outside the sample, in the upper part on the map. Even 
though different in details, the gross features of the density 
field are remarkably similar to those obtained by Zaroubi, 
Hoffman & Dekel (1999) from the application of the same 
formalism to the Mark III catalog. This is an outstanding 
result considering the different ways the two catalogs were 
constructed and the pecuiar velocities measured. 

Fig. 6 compares a higher resolution map of the density 
field recovered from the ENEAR data (left panel) to the den- 
sity field reconstructed from the PSCz redshift catalog (right 
panel; Branchini et al. 1999). Both maps are along the Su- 
pergalactic plane and were reconstructed using a 900 kms -1 
smoothing radius. The shaded area in the left panel indicates 
regions where the error is less than 0.45. Even though differ- 
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Figure 7. The G12 reconstructed velocity field in the Supcrgalac- 
tic plane is displayed as flow lines that start at random points, 
continue tangent to the local velocity field, and are of length pro- 
portional to the magnitude of the velocity at the starting point. 



ent in detail the similarities between the density fields are 
striking lending further credence to the reality of the struc- 
tures observed in the mass distribution. Note that with the 
higher resolution some structures become resolved. One can 
clearly see the Local supercluster at the center of the map 
and that both the GA and PP may consist of different sub- 
structures. 

The velocity field along the Supergalactic plane is pre- 
sented in Fig. 7, showing the existence of two convergence 
regions which roughly coincide with the locations of the GA 
and PP. 



4.3 Bulk Velocity 

The velocity field has been fitted using a monopole, dipole 
[i.e. bulk flow) and quadrupole {i.e. shear) expansion within 
spheres of radii ranging from 1000 to 6000 km s _1 . The three 
Cartesian components of the bulk velocity (in Supergalactic 
coordinates) and its absolute value (Vb) are shown in Fig. 8 
as a function of the depth over which the fit has been done. 
The plots present the bulk velocity of the WF field and of 
an ensemble of 10 CRs. The plot of the absolute value of the 
bulk velocity contains also the mean and standard deviation 
calculated over the ensemble of the CRs. Note that the mean 
Vb of the CRs is higher than its WF value. This result is 
expected as the WF attenuates the velocity field with the 
depth, as the observational errors become more dominant. 

The amplitude of the bulk flow measured from the 
reconstructed three-dimensional velocity field ranges from 
V B = 300 ± 70kms _1 for a sphere of R = 20/i~ 1 Mpc to 
160 ± eOkms -1 for R = 60 /i _1 Mpc. This value is in good 
agreement with that obtained from a direct fit to the ra- 



Figure 8. The bulk velocity fit of the reconstructed velocity field 
as a function of depth: The solid line corresponds to the WF 
field and the dashed lines correspond to an ensemble of 10 CRs. 
The four panels show the (Super Galactic) x,y and z components 
and the amplitude of the bulk velocity. The bottom-right panel 
presents also the mean amplitude taken over the CRs and the 
error-bars are the standard deviation around it. Note that this 
mean value is expected to be larger than the amplitude of the 
WF bulk velocity. 



dial peculiar velocities (da Costa et al. 2000b). This result 
disagrees with the bulk flow determination from the Mark 
III survey where the amplitude is roughly twice that of EN- 
EAR with a comparable scatter (Zaroubi et al. 1999) but 
comparable to that measured from the SFI sample. 



4.4 Large Scale Tidal field 

An alternative way of describing the velocity field is to de- 
compose it into two components, one which is induced by 
the local mass distribution and a tidal component due to 
mass fluctuations external to the volume considered. Here 
we follow the procedure suggested by Hoffman (1998a,b) and 
more recently by Hoffman et al. (2000). The key idea is to 
solve for the particular solution of the Poisson equation with 
respect to the WF density field within a given region and 
zero padding outside. This yields the velocity field induced 
locally, hereafter the divergent field. The tidal field is then 
obtained by subtracting the divergent field from the full ve- 
locity field. Fig. 9 shows the results of this decomposition 
applied to the ENEAR survey, where the local volume is 
a sphere of 80/i _1 Mpc centered on the Local Group. The 
plots show the full velocity field (upper left panel), the di- 
vergent (upper right panel) and the tidal (lower left panel) 
components. To further understand the nature of the tidal 
field its bulk velocity component has been subtracted and 
the residual is shown in the lower right panel. This residual 
is clearly dominated by a quadrupole component. In princi- 
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Figure 9. Tidal field decomposition of The G5 reconstructed 
velocity field in the Supergalactic plane is displayed as flow lines. 
The top-left panel shows to the full velocity field. 

pie, the analysis of this residual field can shed light on the 
exterior mass distribution. 

For the ENEAR catalog we find that the local dynamics 
is hardly affected by structure on scales larger than its depth. 
For this sample not only the bulk velocity at large radii is 
small but so is the rms value of the tidal field estimated to 
be of the order of 60kms~ 1 . This is in marked contrast to 
the the results obtaine from the analysis of the Mark III sur- 
vey which yields a much stronger tidal field, pointing (in the 
sense of its quadrupole moment) towards the Shapley con- 
centration. For Mark III the tidal field contributes roughly 
third of the total bulk velocity (~ 200km/s). 



5 CONCLUSION 

In the first part of this paper the maximum-likelihood 
method (Zaroubi et al. 1997) has been used to measure the 
mass-density power spectrum from the newly completed EN- 
EAR early-type redshift-distance survey. The method as- 
sumes that the galaxy peculiar velocities satisfy Gaussian 
random statistics and that they are linearly related to the 
mass-density field. The initial fluctuation power spectrum is 
assumed to be CDM-like, with or without COBE normal- 
izations. In addition the measured peculiar velocities error 
are assumed to be proportional to the distance with some 
thermal component to account for the nonlinear evolution of 
high-density environment in which the early-type galaxies 
reside. 

General results that are valid for all the models used in 
the analysis, and are independent of the detailed parameter- 



ization and normalization used in each model, can be sum- 
marized as follows. The amplitude of the power spectrum at 
k = O.lhMpc- 1 is P(k)ft 1/2 = (6.5 ± 3) x 10 3 ( ftMpc" 1 ) 3 
yielding r/g = l.llg'35. For the family of COBE-normalized 
CDM models the following range of parameters was con- 
sidered: fl < 1; 0.4 < h < 1; and n = 1. Within this 
range we have obtained a constraint on a combination of 
the parameters fl and h which can be approximated by fl « 
(0.38±0.08)ft _1 ' 3 for ACDM, and fl « (0.52 ± 0.083)/i~ SS 
for OCDM. For h = 0.65, ACDM yields Q = 0.5- 0.8. Simi- 
lar constraints are obtained from the analysis of the generic 
T-models, independent of the COBE normalization. We find 
that the power spectrum amplitude and shape parameter are 
constrained to be r/g — 1.0l ' 2 g an d T > 0.18, with larger 
values of T (> 0.4) being more probable. We point out that 
these constraints are consistent with the results obtained 
from a similar analysis of the Mark III and the SFI pecu- 
liar velocity catalogs. This agreement is encouraging since 
it shows that the results are robust and independent of the 
sample used. 

Examination of the \ 2 /d.o.f. for the most likely COBE- 
normalized models shows that their values are of the order 
of 0.93. These values are about 2a away from the preferred 
value of 1. However, this should not be too alarming as many 
of the models within the errors have x 2 /d-o.f ~ 1. The 
X 2 /d.o.f for the best-fit T-model is 0.99. 

As pointed out by previous papers that have analyzed 
the PS derived from peculiar velocity data (Zaroubi et al. 
1997, Freudling et al. 1999), the constraints on rjg and T arc 
considerably higher than those obtained from other types 
of analyzes including peculiar velocity data (Borgani et al. 
1997, 2000), cluster abundances, and the galaxy power- 
spectrum (Efstathiou et al. 1992; Sutherland et al. 1999). 
They are also not consitent with those obtained by com- 
bining the results from high redshift supernovae type la 
(Perlmutter et al. 1999) and the CMB data (Efstathiou 
et al. 1999) which yields values of fl « 0.25 ± 0.15 and 
A « 0.65 ±0.2. Furthermore, assuming a linear galaxy-mass 
relation the value of rjg obtained from the present analysis 
would imply /3 = TO or a /3i ~ 1.4 (e.g.Willmer, da Costa 
& Pellegrini 1999; Sutherland et al. 1999), where the sub- 
script refers to IRAS galaxies, at least a factor of 2 larger 
than those derived from a velocity-velocity comparison of 
the IRAS 1.2 Jy gravity field and the Mark III (Davis et al. 
1996), SFI ( da Costa et al. 1998) and ENEAR (Nusser et 
al. 2000) all leading to (3i ~ 0.5. These values are also con- 
sistent with those derived from small-scale velocities (Fisher 
et al. 1995). 

There are many possible explanations for the above dis- 
crepancies. One possibility is that all other analyses have 
somehow conspired to produce consistent results but yet 
incorrect interpretation. Even though at first glance this 
seems unlikely, the recent results from the CMB ballon 
experiments Boomerang (Bernardis et al. 2000; Lange et 
al. 2000) and MAXIMA (Hanany et al. 2000; Balbi et al. 
2000) have shown that the height of the second peak in the 
CMB angular power-spectrum is consistent with higher val- 
ues of ft. From their most likely models these authors derive 
fl = 0.4 - 0.8 

It is important to point out that the method is very 
sensitive to the assumed error model which can add or su- 
press power. It also implicitly gives a high weight to nearby 
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galaxies, likely to be slow rotators or low velocity dispersion 
systems, for which the measurements and the distance rela- 
tions are the least reliable. However, tests have shown that 
these effects are not very important for the present data set. 
Another potential problem arises due to the rapid decrease 
of the weight with distance, the effective volume of the cur- 
rently available catalogs is small and the shape of the power 
spectrum is poorly constrained, as illustrated by the case of 
the F-model. All these factors may impact on the reliability 
of the constraints obtained from the PS analysis. 

Finally, one or more of the theoretical model ingre- 
dients could be inaccurate, e.g., power spectrum assumed 
shapes, Gaussianity of the distribution; or even some inher- 
ent bias in the method itself that has eluded the extensive 
numerical tests carried out with the data and mock sam- 
ples (e.g.Freduling et al. 1999). In fact, through an eigen- 
mode expansion of the Mark HI and SFI galaxy catalogs, 
Hoffman and Zaroubi (2000) have conducted a mode-by- 
mode goodness-of-fit analysis. They found that when the 
surveys are analyzed with their corresponding CDM most 
likely models, there is a systematic inconsistency between 
the data and the 'best-fit' models suggesting either a generic 
problem in the peculiar velocity data sets or the inadequacy 
of the theoretical or error models. Unfortunately, however, 
the analysis has not been able to point to the exact source 
of inconsistency. 

Finally, in this study we have also performed, given the 
most probable power spectrum, a Wiener reconstruction of 
the density and velocity fields. The maps shown here have 
1200kms _1 and 900kms _1 Gaussian resolution and they 
are limited to the Supergalactic plane. The main features 
shown are similar to the features in the IRAS reconstruc- 
tion, corrected for perculiar velocities. The constrained real- 
izations allow us to estimate the point-by-point uncertainties 
in the recovered maps. In terms of their recovered density 
fields ENEAR, SFI and Mark III mostly agree. However, 
they do differ in the velocity fields. ENEAR shows no signif- 
icant tidal component which contributes about half of the 
Mark III local bulk velocity. This tidal field accounts for 
the very different bulk velocities obtained from ENEAR and 
Mark III, with SFI situated in between these surveys. The 
results suggest that volumes of 60 — 80/i _1 Mpc are essen- 
tially at rest relative to the CMB and that the Local Group 
motion is primarily due to mass fluctuations within the vol- 
ume sampled by the existing catalogs of peculiar velocity 
data. 
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